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New Discretization of Complex Analysis: 
The Euclidean and Hyperbolic Planes 

Abstract. Discretization of Complex Analysis on the Plane based on the 
standard square lattice was started in 1940s. It was developed by many people 
and also extended to the surfaces subdivided by the squares. In our opinion, 
this standard discretization does not preserve well-known remarkable features 
of the Completely Integrable System. These features certainly characterize 
the standard Cauchy Continuous Complex Analysis. They played a key role 
in the great success of Complex Analysis in Mathematics and Applications. 
Few years ago we developed jointly with I.Dynnikov a New Discretization of 
Complex Analysis (DCA) based on the two-dimensional manifolds with col- 
ored black/white triangulation (see J^). Especially deep results were obtained 
for the Euclidean plane with equilateral triangle lattice in the works J2l [7^ . 
Our approach preserves a lot of features of Completely Integrable Systems. 
In the present work we develop a DCA theory for the analogs of equilateral 
triangle lattice in the Hyperbolic plane. This case is much more difficult 
than Euclidean. Many problems (easily solved for the Euclidean Plane) are 
not solved here yet. Some specific very interesting "dynamical phenomena" 
appear in this case: for example, description of boundaries of the most fun- 
damental geometric objects (like the round ball) leads to dynamical problems. 
Mike Boyle from the University of Maryland helped me to use here the meth- 
ods of symbolic dynamics. 



Introduction. History. We do not discuss here "geometric" discretiza- 
tions of conformal mappings started in early XX century-see survey and his- 
tory in [S]. By the way, specific topological properties of surfaces subdivided 
by the squares, was started by the combinatorial geometry/topology group 



^Sergey P. Novikov, IPST and MATH Department, University of Maryland, Col- 
lege Park MD, USA and Russian Academy of Sciences (Steklov Math Institute and 
Landau Institute for Theoretical Physics, Moscow), e-mails novikov@ipst.umd.edu, 
snovikov@mi.ras.ru; This work is partially supported by the Russian Grant in the Nonlin- 
ear Dynamics. 



1 



in the Steklov Institute in 1980s by the suggestion of the present author (see 
[B]). This group called this object "The Quadrillage" . I formulated these 
problems under the influence of the Statistical Physics of Lattice Models, 
after discussions with A.Polyakov. It is mentioned also in my Encyclopedia 
volume number 13, Topology- 1. -see [7]. In the works of the group [5] made in 
last decade such surfaces subdivided into squares are called " quad-graphs" . 

Our goal is to discretize Cauchy-Riemann operator 9 as a Linear 
Difference First Order Operator on Triangulated Surfaces. Let us 
say at the very beginning of this article that Neither in the Standard 
Approach nor in our New Discretization The Discrete Analogs of 
Holomorphic Functions form a Commutative Associative Ring (at 
least naturally). 

The Discretization of Complex Analysis as a Linear Difference Operator 
was done for the square lattice in by Lelong-Ferrand. Her work was 
published in 1944 (see [2]). In the work [3] Duffin developed or rediscovered 
this approach. He also extended it to the rombic lattices later. 

Duffin claimed that another mathematician already discovered this idea 
before and published his work in 1941 in the notes of some small provincial 
South American University. Our attempts to find this obscure edition failed, 
so we continue to quote Lelong-Ferrand as a first author who invented this 
idea. 

In the standard approach Discrete Analog of the Cauchy-Riemann Oper- 
ator d acts on the C-valued functions ip of vertices in the square lattice on 
the plane: 

Qsquarei>{m, u) = il){m, u) + z?/;(m + 1, n) - ■il){m + 1, n + 1) - iil){m, n + l) 

By definition, every solution to the equation Qsquarei^ = is a d-holomorphic 
function. 

A lot of people work here now developing this approach (see, for example, 
in[5]). 

Let us present a hint of the idea: Why we are not satisfied by the standard 
discretization? Of course, every continuous model admits many difference 
approximations which converge to it in the continuum limit. 

For the continuous systems with huge hidden algebraic symmetry we wish 
to find optimal discretization preserving as much symmetry as possible. As we 
already know, there exists no natural discretization preserving multiplication 
of holomorphic functions. All attempts to invent multiplicative structure with 
good properties failed. It became clear many years ago. 
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Let us point out following: 

1. In the standard approach discrete analog of the Cauchy-Riemann 
operator is in fact a second order difference operator: two lengths are involved 
in the sum (the lengths of sides and diagonals). 

We wish to have a first order difference operator as a natural 
discrete analog of d. Is it possible? The answer is YES for the 
Equilateral Triangle Lattice. Corresponding theory of the differ- 
ence first order "Triangle Operators" on the simplicial complexes is 
developed in the Part 1 in more general form. We define also non- 
standard discrete analogs of GL„-Connections for the triangulated 
n-manifolds. 

In the Part 2 we develop a general theory of discrete analogs of holomor- 
phic functions (the d-holomorphic functions) on the 2-manifolds with dis- 
crete analog of conformal structure (DCS). We define DCS as a colored 
black/ white triangulation. 

The d-analogs of Liouville Principle and Maximaum Principle will be 
clarified here for the general 2-surfaces with DCS. 

Part 3 is dedicated to the theory of d-holomorphic functions on the Equi- 
lateral Triangle Lattice. 

2. Continuous 2D Laplace operator admits a natural factorization 

A = aa 

Nothing like that exists on the square lattice. 

However, Equilateral Triangle Lattice is much better for this-see below. 
Problems: 

1. How to construct Discrete Analogs of Holomorphic Polynomials without 
multiplication? 

People did it with great efforts after many years in the standard approach. 
Their constructions are non- canonical and non-unique. 

It is very easy in our new approach. Our construction is obvious 
and canonical. 

2. How to construct rational functions without multiplication? 

In particular, in the standard continuous complex analysis rational func- 
tions are especially important because the function 1/ z (the Cauchy Kernel) 
is also a unique fundamental solution to the equation 

d{%l){z)) = 2m6{z) 
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decreasing at infinity? 

The answer to this question is also positive and canonical in our 
approach. It was missed in our first work with Dynnikov [9] and 
made by Grinevich and R.Novikov later in |10| . 

Part 4 is dedicated to the Hyperbolic (Lobachevski) Plane H. 

Equilateral Triangle Lattice in H is defined simply as a triangulation of 
plane such that p triangles hit every vertex and p > 7. We need p = 2s. For 
all s > 4 we have a triangulation of Lobachevski Plane H such that it admits 
a black/white coloring. Let us concentrate our attention on the case s = 4. 

At the moment we can construct neither d-analog of polynomials nor 
d-rational functions. 

Even boundary of round ball Dr of radius r E Z is a. complicated object. 
How to describe it? How many points does it contain? We succeeded in the 
solution of this problem using methods of symbolic dynamics. Mike Boyle 
from the University of Maryland helped us. The boundary dDr contains 
approximately A'' points where A = 2 + ^3- 

It turns out that for the reconstruction of d-holomorphic function in the 
ball Dr+i we need to know values of the boundary function in [n(r + 1)/2] + 1 
points. Here n(r) is a number of vertices in the boundary. It exactly corre- 
sponds to the continuous case where only set of coefficients at the nonneg- 
ative powers of z (i.e. one half plus one Fourier coefficients of function on 
the boundary circle) determines completely our function inside. The main 
unsolved problem is: 

Problem. How to find a bounded basis of d-holomorphic func- 
tions on Hyperbolic Plane? 

Part I. Definitions. Discrete GLn-Connections. B/W manifolds. 
The First Order Triangle Operators. 

Take a simplicial complex M. We assume that this complex is given with 
canonical metric where every simplex is a standard unit linear subsimplex in 
euclidean space with natural euclidean metric. We fix a family of ra-simplices 
such that every vertex belongs at least to one simplex of that family. Fix 
also collection of nonzero numbers Btip 7^ for all T G X and P E T. 

Definition 1 We call following operator the Triangle Operator associ- 
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ated with family X: 

PeT 

It maps functions of vertices into functions of simplices T G X. Our co- 
efficients normally belong to R or C. The operators are the first order 
linear difference operators. 

Now we define a Discrete Version of GL„-Connections. 

Definition 2 Let X be a family of all n-simplices in n-manifold M. We call 
the equation Qip = Discrete Differentially- Geometrical GLn- Connection. 
The coefficients hx-.p are defined in every n-simplex up to nonzero factor. So 
every DC- Connection is defined by the set of ratios 

Tiiis discretization is different from tlie standard Wilson Discretization 
used by pliysicists studying Yang-Mills fields; There is no natural way to 
select compact holonomy groups in our approach. 

The theory of discrete GLn connections was constructed in the works [H 
H] , some first ideas appeared in [H [9] in connection with completely integrable 
systems. We are working only with scalar-valued functions. At the same time 
we are going to define a Nonabelian Curvature and Holonomy with values in 
the group GL^ for the n-manifolds. 

As we demonstrated in the [Ij, two different Holonomy Representations 
can be constructed here: Abelian ("Framed Abelian") and Nonabelian. We 
need here only Nonabelian Holonomy. 

What is a Nonabelian Discrete Curvature? A Nonabelian Holonomy 
group is defined along the Thick Paths. 

Definition 3 We call by Thick Path any sequence of n-simplices < TiT2...Tk > 
such that intersection Tjf]Tj^i is exactly n — 1- dimensional face Aj, and 
Aj 7^ Aj^i. We say that Thick Path is closed if T^ = Ti. Its length is equal 
tok-l. 

Therefore in every simplex Tj of the thick path exactly two different n — 1- 
dimensional faces are fixed: Aj C T^+i is called an "In-face" in Tj+i, and 
A J C Tj is called an "Out-face" in Tj. 
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Let us define a Parallel Transport of the n- vector ip consisting of values 
in all vertices of the in-face Ai C Ti, along the thick path < Ti...Tk. Solving 
equation Q = 0, we define this function in all vertices of Ti. So we know 
values of in all vertices of the out-face in Ti equal to the in-face in T2. After 
that we are doing the same procedure for T2 and so on. Finally we define a 
"Transported value" of n- vector ip in the out-face of the simplex T^. This map 
is linear. For the closed paths we are coming to the Nonabelian Holonomy 
Representation of the Semigroup of the closed thick paths into the 
Group GLn. 

Definition 4 We call by the Nonabelian Discrete Curvature in the vertex P 
a Holonomy Linear Map along the Thick Closed Paths in the Simplicial Star 
St{P). It is enough to know holonomy for all thick closed paths around all 
n — 2=simplices-see Fig 1 for n = 2. 




However, the natural gauge group 

^^fij,Q^g{T)Qf-\P) 

is abelian, / 7^ 0, 7^ 0. There are very interesting specific features here 
realizing non-abelian GLn Connections using only the spaces of scalar func- 
tions of vertices. Some sort of mixing abelian and non-abelian properties 
appears. Details of classification can be found in the work [1]. No theory of 
Characteristic Classes is constructed yet. 

Definition 5 We call special discrete GLn- Connection with bx-.p = 1, a 
Canonical Connection 
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Consider now any n-manifold. Is it possible to color its n-simplices by 
the black and white colors such that for every pair of n-simplices attached 
to each other along n — 1-face, the colors are opposite? 

Definition 6 We call such colored manifolds a B/W manifolds. 

Such coloring exists if and only if every closed thick path consists of even 
number of n-simplices. More general, fix any n — 1-simplex A. We have a 
natural homomorphism 

^2:^mck{M,^)^ Z2 = Z/2Z 

defined as a number of n-simplices in any thick path modulo 2. This map 
also looks like simplest Connection whose Holonomy Group is Z2, but we do 
not realize it by the linear operator. 

Definition 7 Let the map $2 is trivial for all thick closed paths in the star 
of every vertex. We call such triangulation "Locally B/W" . 

In this case following simple lemmas are true: 

Lemma 1 Under this assumptions the map $2 defines a map 02 : 7ri(M) — >■ 

Lemma 2 Curvature of the Canonical Connection is trivial for every locally 
B/W manifold. 

Proof. Fix any n — 2-simplex S in the star St{P). Let exactly s{S) n- 
simplices contain S. The Nonabelian Holonomy with initial values ao, ...,a„ 
in the vertices of n-simplex T„ =< S,Pn-i,Pn > containing S (and around 
it), is determined by permutation A^^^\ Here A{ao, an-2, ctn-i, o.n) = 
(ao, an_2. On, a„_i). So = 1. We see that Nonabelian Holonomy (Cur- 
vature) is equal to 1 if and only if s{S) is even. Lemma is proved. 

Every locally flat Connection (i.e. Nonabelian Curvature is trivial for 
every vertex) defines a homomorphism of topological holonomy vri(M) — >■ 

GLn- 

Lemma 3 Every flat Canonical Connection defines a Topological Holonomy 
Homomorphism. Its image belongs to the subgroup Sn+i C GLn. 
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Proof. The imbedding is following: Take arbitrary distinct n+1 numbers 
Oq, such that Yl^j ~ 0- ^^^^ °f numbers in the vertices of 

initial n-simplex Ti. For every thick closed path Ti,....,Tfc the values of 
Parallel Transport along this path in every simplex Tj consists of the very 
same numbers because Yl^j ~ 0' Connection is Canonical. In the 
final simplex = Ti we obtain permutation belonging to Sn+i which does 
not depend on the choice of initial numbers aj. Lemma is proved. 

Now we are going to consider only Orientable Manifolds M. 

Lemma 4 Let Curvature and Topological Holonomy are trivial for the Canon- 
ical Connection (i.e. it is flat). Then B/W coloring exists globally. 

Proof. Let a closed thick path < Ti, > defines a trivial permutation 
corresponding to the Canonical Connection. We claim that this path consists 
of the even number of n-simplices. Indeed, it follows from the fact that our 
Nonabelian Holonomy is equal to the product of elementary transpositions. 
Every one of them changes orientation. The whole path preserves orientation. 
So our lemma follows. 

Lemma 5 For every manifold M with flat Canonical Connection the equa- 
tion Qtp = defines n-dimensional linear space of Covariant Constants. 

Proof obviously follows from the fact that every initial values in the n — 1- 
simplex Delta defines a covariant constant correctly because our connection 
is flat. 

So we fix a manifold M with B/W structure. We always assume that 
our Canonical Connection is flat. It is always true for simply connected 
manifolds. It is locally fiat because B/W structure is given. For non-simply- 
connected case B/W structure implies that Nonabelian Holonomy maps 7ri(M) 
into Sn+i- So Canonical Connection became globally fiat on some finite cov- 
ering with number of sheets no more than n\. We define d- analog of holo- 
morphic functions in every B/W manifolds (see below). 

How to construct Continuous Limit for our Discretization of 
Complex Analysis? 

For that we are doing the following procedure: Let n = 2. Construct 
a special Covariant Constant /o whose values in every triangle are 1, (, 
where = 1. It is unique up to the group 5*3 permuting vertices of the initial 
triangle. Apply following gauge transformation to the Canonical Connection: 
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We still have two-dimensional space of covariant constants for this connection 
but in this gauge one of covariant constants became an ordinary constant. 
We work over the extended field containing (. For the case R and extended 
field C we used this construction to get Continuum Limit of our theory to 
the ordinary complex analysis. For M = B? exactly one half of our theory 
converges to the ordinary complex analysis. The second half diverges for 
small triangulations. 

Now we start to construct New Discretization of Complex Anal- 
ysis. 

Let M be any B/W-manifold. Following two families play fundamental 
role in our theory: Xi = b (all black simplices) and X2 = w (all white 
simplices)-see Fig 2. 




Fig 2 K=0 



Definition 8 The operators Q'', are called Black and White Triangle Op- 
erators. The union of this families hl^w is complete. Using the Canonical 
Connection, we define an operator Q = @ . 

Following lemma is very simple: 

Lemma 6 The operators Q^, Q*", Q have following properties: 

L = Q*Q = 2Q^*Q^ = 2Q'"*Q'" 
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Let A = dd* he an ordinary scalar simplicial Laplace- Beltrami Operator, nip 
is the number of edges entering vertex P. For 2-manifolds we have Q*Q = 
-2A + Snip. 

Proof follows from thew elementary calculation. 

Definition 9 Consider any function ip such that Q^ip = 0. Forn = 2 we call 
such functions (i-holomorphic. We call solutions to the equation Q'^ip = 
d-antiholomorphic. 

Now we assume that the Canonical Connection is flat. For every black 
triangle T we define a unique covariant constant EipiT) whose values in T 
coincide with -0. 

Definition 10 We call Eip an Evaluation of function ip in the black 
triangle T. 

As a conclusion, we can view every d-holomorphic function ei- 
ther as a real-valued function ip on the vertices or as a i?^-valued 
function Eip(T) of the black triangle T. 



Part II. D-holomorphic functions on 2-manifolds. Liouville Prin- 
ciple and Maximum Principle 



As above, we consider triangulated orientable B/W manifolds M with 
fiat Canonical Connection. Our field here is R. Following general properties 
of d-holomorphic functions were found in [1] : 

Theorem 1 (The Liouville Principle). For the compact closed 2-manifold 
every d-holomorphic function is a covariant constant. 

Proof is based on the "instanton phenomenon" for the quadratic func- 
tional {Lif)., Ip) if L is factorizable L = Q*Q. Let the global minima ip in 
the Hilbert space L2 is a zero mode Lip = 0. The order of Euler-Lagrange 
equation drops twice reducing to the "Self-Duality Equation" which is in 
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our case simply the d-analog of the Cauchy-Riemann Equation. Let us re- 
mind here that This Property is the most fundamental feature of 
the standard continuous Complex Analysis unifying it with The 
Completely Integrable Systems. The proof of our theorem is following: 
= implies {Q^)*Q''ip = implies = Q*Qijj = implies {Lijj, tjj) = 
implies {Qip, Qip) = implies Qip = 0, so our theorem is proved. 

In principle, the famous Instanton Phenomena discovered by Polyakov, 
Belavin, Schwarz and Tiupkin in 1974 in the Nonlinear Variational Calculus 
(Yang-Mills Field Theory), is based on the similar arguments. 

Let D be a bounded domain in M consisting of black triangles. 

Definition 11 We call T & D a boundary triangle if some of its vertices 
belong to black triangle not belonging to the domain D. 

Theorem 2 (The Maximum Principle) For every d-holomorphic function ip 
in D the set of covariant constants Eip{T) G for all T E D is contained 
in the convex hull of the image of boundary triangles. 

Proof can be found in [1]. 

In the next chapters we consider specific results obtained for the Euclidean 
(Part 111) and Lobachevski or Hyperbolic Planes (Part IV). 

Part III. Equilateral Triangle Lattice in R"^. 

The standard regular equilateral lattice in R^ (see Fig 3) admits two basic 
(unitary) shift operators ti,t2 acting on vertices 

ti(m, n) = (m + 1, n),t2{m, n) = (m, n + 1) 
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All six elementary shifts , (titg ))^^ have equal length. The operators 
Q'^, are equal to 

= 1 + tl + t2, = (Q')* = 1 + tr' + ^2 ' 

They map the space of functions of vertices into itself and are adjoint to 
each other. Our field here is R. In the works [H [9] we defined first time 
black and white triangle operators and invented the idea of new type discrete 
GL2 connection. We proved also that every difference second order self- 
adjoint operator L = a + bti + ct2 + (itit^^ + (adjoint) admits Laplace- 
type factorization L = Q*Q + W where Q is some triangle operator Q = 
u + vti + wt2. and W, a, b, c, d, u, v, w are some real functions on this lattice. 
In particular, —A + 9 = Q^Q^ for the standard Laplace-Beltrami operator. 
A theory of discrete completely integrable systems based on the discretized 
second order operators, was started. 

Definition 12 We call d-holomorphic function ip polynomial (i.e. ip E Polk) 
of degree k if Q''ip = and {Q'")''+^ip = 0. 

As it was established in [1], these functions have a /c-polynomial growth 
at infinity in R^. They are completely determined by their values in any 
standard A;-triangles , black from inside, with 2k + 2 points in each edge 
(see Fig 4). 
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So the dimension of the space Polk of /c-polynomials is 2k + 2. We can 
choose the canonical /c-polynomials 4'Tk,a, a = 1, 2, 3. They are equal to zero 
in Tk everywhere except one boundary edge a where they have values ±1 . 
Sum of them belongs to the space Polk^i. 

^ i)k,a e Polk^i 

a 

Therefore for the choice of basis we have to select pair of edges in Tk- 

Following theorem ( [1] ) give natural analog of polynomial approximation 
of holomorphic functions leading to the " d- Taylor Series" : 

Theorem 3 For every function ip , d-holomorphic in , and every canonical 
triangle Tk, there exists a unique k-polynomial (j) such that ip — cf) is identically 
equal to zero in the triangle Tk- 

For the Taylor decomposition of we have to choose an increasing se- 
quence of triangles T^, k = 0, 1, 2, ... and pair of boundary edges in each of 
them. The choice of such sequence is non=canonical. 

How to get analog of the Cauchy formula? 

We need to construct a " Cauchy Kernel" (d-analog of 1/z) satisfying to 
equation: 

Q''G{x,y) = 5{x-y) 

where x,y are points in the lattice x = {m,n),y = {m',n'), and difference 
operator Q'^ acts on the variables x = {m,n). Having any such function, we 
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consider any bounded domain D and (i-holomorphic function i/j in D. Let us 
extend this function to the function ip such that 

ijj{x) = tp{x),x e D; ^{x) = 0,x E R^minusD 

We have 

J2[Q'i^iy)]Gix-y)=^lj{x) 

y 

for all X G D. Let us point out that the function Q^ip{x) is concentrated 
along the " boundary strip" , so it is really analog of the Cauchy formula. In 
the work [1] we constructed following "hyperbolic-like" Cauchy Kernel (see 
Fig 6). We call it Pascal Triangle. It is equal to zero outside of infinite 
triangle and has exponential growth inside of it. 

.0 



1 y///A'////W^ 



/////Ay// 


//)vW/ 


vM*/////y\ 





\///Av 

Fig 5 "Pascal Triangle" G(x) x=(m,n) 

As it was pointed out by P.Grinevich and R.Novikov in [10], much better 
Green function can be obtained simply by the Fourier Transform 

G{x) = (const) / / dkidk2e'"''''e'''''^{l + exp{ifci} + expjz/ta})"^ 
Jo Jo 

where x = (m, n). 

They proved following theorem: 

Theorem 4 This integral converges. This Green Function decays at infin- 
ity as c/d{P) where d{P) is equal to the distance of this point to zero. In 
particular, every d-holomorphic function on the whole equilateral triangle lat- 
tice whose grouth at infinity is no more that polynomial, is a d-holomorphic 
polynomial. 
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Conclusion. This result gives a unique function decaying like 
l/|x| at \x\ oo. All functions (Q'^)''G{x) has a growth like 
Therefore we have also analogs of rational function without any of 
multiplication (which does not exist). 

Part IV. D-holomorphic functions on The Hyperbolic (Lobachevski) 
Plane. 

Every triangulated plane such that more than six triangles (edges) hit 
every vertex, can be viewed as a negative curvature plane. An Equilateral 
Triangle Lattice we get in the case then this number is the same for all vertices 
(i.e. mp > 7). We need also B /W structure so our Equilateral Triangle 
Lattices in with d-conformal structure (i.e. with B/W coloring) 
correspond to m = 8, 10, 12, .... We consider here only the case m = 8. 

Let us point out that the triangles of this lattice cannot be made arbitrary 
small in the standard Lobachevski metric (their size is fixed in by the 
number m) but our domain D can be made arbitrarily large. 

For example, for every finite set of vertices K C and positive integer 
r G we define a domain D^^r consisting of vertices x G D^^r such that 
distance d{x, K) is no more than r. We measure distances between vertices 
(sets) counting minimal number of edges needed for joining them. The sim- 
plest important domains of that kind correspond to the cases: l.i^=vertex 
(we call it "standard ball" D,.); 2.K=triangle T i.e. K=?) vertices of T. We 
denote it Dt^t- It is also like a ball with center in the center of triangle. For 
K = Dk',1 we have DK,r = Dk'^t+i if K is connected (i.e. no jumps of the 
length more than one are needed to reach one point from another). 

There is very big automorphism group mapping this Lattice into itself 
but this group is non-commutative. It does not contain big enough commu- 
tative subgroups, so nothing like Fourier transform exists here. We cannot 
construct good enough Green function. Our operators Q'', map space of 
functions of vertices into the space of functions of black and white triangles 
correspondingly. Therefore we cannot iterate them. So we cannot construct 
analogs of polynomials here. 

Problem: How to construct basis of d-holomorphic functions 
'4)i{x) in Hi which are globally bounded in all space? 

In the continuous case our space is realized as a unit disc : \z\ < 1. 
We have basis of bounded holomorphic functions z^,k = 0,1,2,.... Every 
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rational function with poles outside of unit disc is bounded in it. 

Easy to construct some d-holomorphic functions ) equal to zero 

inside of the r-ball Dr (i.e. for x G Di,l < r) and equal to zero along the 
boundary dDj. except of the specific place P C dDj. looking like in Fig 6, a 
and b: 



1-11 1-1 




Fig 6 P=...wbbw... ...wbw..=P 

Its extension to the external domain is non-unique. Can we construct 
a globally bounded extension?. 

Let us introduce class of right-convex paths. 

Definition 13 We call oriented path consisting of edges Right Convex if it 
bounds two or three triangles only from the right site in every vertex 

We are coding all right convex paths by the words in two symbols w, b 
assigning to every edge letter w or b depending on which color has triangle 
from the right site of this path.-see Fig 7 




bb WW 



Fig 7 c) d) 
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Theorem 5 Let ip be a d-holomorphic function and R is a full set of zeroes 
ip{x) = 0. Consider boundary of its complementary domain, i.e. the set 
Dn^iminusR, which is a set of points-closest neighbors of zeroes. Every con- 
nected component of this boundary set is right convex choosing orientation 
such that corresponding component of zeroes lies inside. 

Proof. In order to have zero value of ip in the point x we should have following 
values in the neighboring point outside of zero set (see Fig 8, a and b). 



C 






Fig 8 ^) ^) 

No other possibilities exist because of the equation Q^ip = 0. It means 
in this case that c 7^ implies either d 7^ or (i = in the Fig 8. In the 
first case our boundary contains only two triangles inside. In the second case 
{d = 0) our boundary contains three triangles inside. Anyway, it is convex. 
Our theorem is proved. 
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A maximal right convex set has a boundary path coded by the infinite 
sequence ...bbbbb.... For every point x and direction-edge / started in x and 
having black triangle from the right, we uniquely construct such maximal 
path ...bbb... and denote it by jx,!- Easy to construct d-holomorphic function 
4'x,i{y) such that ip = Oin the domain to the right of the path 'jx,i, and ip{x) = 
±1 along the path 7a. i-see Fig 9. Its continuation to the complementary 
domain is non-unique. We can define it from the requirement that the growth 
is minimal if distance d{y,'~fx,i) oo, but this definition is noneffective. 
How to find this growth? What is a "minimal" function ijjx,i'^- This 
function and its group shifts give basis in the space of all d-holomorphic 
functions. 

Every right-convex path 7 can be shifted to the left side by the distance 
equal to one. We get a new path T{'-f). Between these two paths we have 
a strip containing all triangles T having at least one common vertex with 
path 7. New path T(7) consists of edges opposite to vertices belonging to 7. 
Only triangles having exactly one common vertex with 7 participate in this 
construction (see Fig 10). 



Fig 10 




Lemma 7 New path T{'~f) is also right convex. Its word can be obtained 
from the word describing 7 by the following procedure: take every pair of 
neighboring letters in 7 and put between them words written below 

bw bwbw, wb wbwb 

bb bwb, WW — ^ wbw 
After that delete old letters. What remains is exactly a word T{'j). 
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Proof easily follows from the picture (see Fig 11). 




a) b) 
bw — bwbw bb — bwb 




This map can be studied by the technic of symbolic dynamics. Mike 
Boyle from the University of Maryland helped me a lot with this business. 
We introduce new letters in order to describe this map. Let bw corresponds to 
the symbol Wb, bb corresponds to bb wb corresponds to by^ and ww corresponds 
to Wui- Written in these symbols, our map T has a form 

Wb WbbwWbWro-,Wyj byjWbWw 

and 

bw byjWbbyjbb-, bb Wbbwbb 

After abelianization and replacing product by sum, we are coming to the" Perron 
matrix" A whose largest eigenvalue is A = 2 + v^- So we proved following 

Theorem 6 The size of right convex path 7 shifted to the left side by the 
distance one, increases asymptotically by the factor A = 2 + ^/S, so we have 
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|T(7)| ~ (2 + -\/3)|7|- In particular, this is true for the boundary of r -ball 
|7| = \dDr\ ~ X'~''\dDk\,k> 1 

Example 1 We have \dDr\ = 8, 32, 120, 448, 1672, ... /or r = 1,2,3,4,5,.., 
so our asymptotic formula is practically exact for r > 3. 

How to calculate dimension of the space of d-holomorphic func- 
tions in the ball How many data on the boundary dDr are 
needed to recover d-holomorphic function in D^? 

Consider strip between Dr and D^+i-see Fig 12, r=l. 




Fig 12 

Every letter b in dDr+i defines black triangle in the strip touching dDr 
in one vertex. Every letter w in dDr defines black triangle in the same strip 
touching dDr+i in one vertex. So total number of equations Q^ip = in 
this strip is equal to Br+i + Wr-, the numbers of black and white letters in 
the boundaries dDr+i, dDr. Unifying all strips with A; < r + 1, we get total 
number of equations Eqr+i : Q'^ip = 0: 

Eqr+i = {Br+l + Wr) + ... + {B^ + Wi) + Bi 

For the number of points in -Dr+i we have 

Nr+l = 1 + Bi + Wi + ...+ Bk + Wk + ... + Br+l + Wr+1 

and \dDk\ = Bk + Wk- So we have Eqr+i = Br+i + Nr — 1. For the number 
of necessary data on the boundary dDr+i we obtain (taking into account 
Wu = Bk) 

Nr+l - Eqr+i = Wr+l + 1 = \dDr+l\/2 + 1 
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Conclusion: Our result that "The Number of Necessary Data 
is equal to [|9i5r+i|/2] + 1" is an exact analog of the Continuous 
case: Only the set of Fourier coefficients on the boundary circle 
corresponding to the exponents exp{m0} with nonnegative n > 0, is 
needed for the reconstruction of holomorphic function in the whole 
disc inside. 
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